Fixed-iV Superconductivity: The Crossover from the Bulk to the Few-Electron Limit 



Fabian Braun and Jan von Delft 
Institut fur Theoretische Festkorperphysik, Universitat Karlsruhe, 76128 Karlsruhe, Germany 

(November 23, 1998) 

We present a truly canonical theory of superconductivity in ultrasmall metallic grains by vari- 
ationally optimizing fixed-iV projected BCS wave-functions, which yields the first full description 
of the entire crossover from the bulk BCS regime (mean level spacing d <C bulk gap A) to the 
"fluctuation-dominated" few-electron regime (d > A). A wave-function analysis shows in detail 
how the BCS limit is recovered for d <C A, and how for d > A pairing correlations become delo- 
calized in energy space. An earlier grand-canonical prediction for an observable parity effect in the 
spectral gaps is found to survive the fixed- A projection. 
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In the early days of BCS theory, its use of essentially 
grand-canonical (g.c.) wave-functions was viewed as one 
of its most innovative, if not perplexing features: the 
variational BCS ansatz for the ground state is a superpo- 
sition of states with different electron numbers, although 
BCS 0] themselves had emphasized that the true ground 
state of an isolated superconductor must be a state of def- 
inite electron number. That this ansatz was nevertheless 
rapidly accepted and tremendously successful had two 
reasons: Firstly, calculational convenience - determining 
the variational parameters is incomparably much simpler 
in a g.c. framework, where the particle number is fixed 
only on the average, than in a canonical one, where a fur- 
ther projection to fixed electron number is required; and 
secondly, it becomes exact in the thermodynamic limit - 
fixed- TV projections yield corrections to the BCS ground 
state energy per electron that are only of order JV -1 , as 
shown e.g. by Anderson Q and Miihlschlegel ||. 

Recently, however, a more detailed examination of the 
range of validity of BCS's g.c. treatment has become nec- 
essary, in the light of measurements by Ralph, Black and 
Tinkham (RBT) |Q of the discrete electronic spectrum of 
an individual ultrasmall superconducting grain: it had a 
charging energy so large [Ec ^ A) that electron number 
fluctuations are strongly suppressed, calling for a canoni- 
cal description, and the number of electrons N within the 
Debye frequency cutoff from the Fermi energy ef was 
only of order 10 2 , hence differences between canonical 
and g.c. treatments might become important. Moreover, 
its mean level spacing d oc iV -1 was comparable to the 
bulk gap A, hence it lies right in the crossover regime be- 
tween the "fluctuation-dominated" few-electron regime 
(d > A) and the bulk BCS regime (d <C A), which could 
not be treated satisfactorily in any of the recent theo- 
retical papers inspired by these experiments: the results 
of including the predictions of parity effects, were 

obtained in a g.c. framework; and Mastellone, Falci and 
Fazio's (MFF) [jl0| fixed- N exact numerical diagonaliza- 
tion study, the first detailed analysis of the f.d. regime, 
was limited to N < 25. 

In this Letter we achieve the first canonical descrip- 



tion of the full crossover. We explicitly project the 
BCS ansatz to fixed N (for N < 600) before varia- 
tionally optimizing it, adapting an approach developed 
by Dietrich, Mang and Pradal [[ll| for shell-model nu- 
clei with pairing interactions to the case of ultrasmall 
grains. This projected BCS (PBCS) approach enables us 
(i) to significantly improve previous g.c. upper bounds 
on ground state energies ; (ii) to check that a previ- 
ous grand-canonical prediction p] for an observable par- 
ity effect in the spectral gaps survives the fixed-iV pro- 
jection; (iii) to find in the crossover regime a remnant 
of the "break-down of superconductivity" found in g.c. 
studies, at which the condensation energy changes from 
being extensive to practically intensive; and (iv) to study 
this change by an explicit wave-function analysis, which 
shows in detail how the BCS limit is recovered for d <C A, 
and how for d 3> A pairing correlations become delocal- 
ized in energy space. 

The model. — We model the superconducting grain 
by a reduced BCS-Hamiltonian which has been used be- 
fore to describe small superconducting grains J@-§| (it 
was phenomenologically successful for d < A [|7flg], but 
probably is unrealistically simple for d > A, for which it 
should rather be viewed as toy model): 
N—l JV-1 
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The ct ± create electrons in free time-reversed single- 
particle- in-a-box states \j, ±), with discrete, uniformly 
spaced, degenerate eigenenergies Sj = jd + e . The in- 
teraction scatters only time-reversed pairs of electrons 
within lud of ef- Its dimensionless strength A is related 
to the two material parameters A and ojd via the bulk 
gap equation sinhl/A = qjd/A. We chose A = 0.22, 
close to that of Al ||. The level spacing d determines 
the number N = 2uj]j/d of levels, taken symmetrically 
around ef, within the cutoff; electrons outside the cutoff 
remain unaffected by the interaction and are thus ne- 
glected throughout. 

Projected Variational Method. — We construct varia- 
tional ground states for H by projecting BCS-type wave- 
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functions onto a fixed electron number N = 2no + 6 [|ll| , 
where no and b are the number of electron pairs and 
unpaired electrons within the cutoff, respectively. Con- 
sidering 6 = first, we take 
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where |Vac) is the vacuum state. Both Vj, the amplitude 
to find a pair of electrons in the levels \j, ±), and Uj, the 
amplitude for the level's being empty, can be chosen real 
p| and obey uj + vj — 1 . The integral over <fi performs 
the projection onto the fixed electron pair number no, 
and C is a normalization constant ensuring (0|0) = 1. 

Doing the integral analytically yields a sum over ( j^ ) 
terms (all products in (|J) that contain exactly no factors 
of VjCj + Cj_ ), which is forbiddingly unhandy for any rea- 
sonable no . Therefore we follow Ref . jll| and evaluate all 
integrals numerically instead. Introducing the following 
shorthand for a general projection integral, 
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Minimization with respect to the variational parameters 
Vj leads to a set of 2no coupled equations, 



Aj)u jVj = Aj( Uj -v ), 
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where the quantities ij, Aj and Aj are defined by 
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We obtain an upper bound on the ground state energy 
and a set of u, 's, i.e. an approximate wave function, by 
solving these equations numerically. To this end we use a 
formula of Ma and Rasmussen JI!| to express any R^'" 3N 
in terms of Ro and all Rq's, and evaluate the latter inte- 
grals using fast Fourier transform routines. 

Next consider states with b unpaired electrons, e.g. 
states with odd number parity or excited states: Un- 
paired electrons are "inert" because the particular form 
of the interaction involves only electron pairs. Thus the 



Hilbert subspace with b specific levels occupied by un- 
paired electrons, i.e. levels "blocked" to pair scattering 
p3|,[7t, is closed under the action of H, allowing us to cal- 
culate the energy, say £b, of its ground state |6) by the 
variational method, too. To minimize the kinetic energy 
of the unpaired electrons in 6} we choose the b singly oc- 
cupied levels, j € B, to be those closest to the Fermi sur- 
face Our variational ansatz for |6) then differs from 

|0) only in that flf is replaced by (HjeB c ]+) I\j? B - 
Thus in all products and sums over j above, the blocked 
levels are excluded (the uj and Vj are not defined for 
j £ B) and the total energy £b has an extra kinetic term 

In the limit d — » at fixed nod, the PBCS theory re- 
duces to the g.c. BCS theory of Ref. @| (proving that the 
latter's iV- fluctuations become negligible in this limit): 
The projection integrals can then be approximated by 
their saddle point values (Tl); since <fi = at the saddle, 
the R's used here are all equal, thus Aj vanishes, the vari- 
ational equations decouple and reduce to the BCS gap 
equation, and the saddle point condition fixes the mean 
number of electrons to be 2no + b. To check the oppo- 
site limit of d ^> A where no becomes small, i.e. the f.d. 
regime, we compared our PBCS results for £q and £\ with 
MFF's exact results |l(J, finding agreement to within 1% 
for no < 12. This shows that "superconducting fluctu- 
ations" (as pairing correlations are traditionally called 
when, as in this regime, the g.c. pairing parameter van- 
ishes ||) are treated adequately in the PBCS approach. 
Because it works so well for d <C A and d ^> A, it seems 
reasonable to trust it in the crossover regime d ~ A too, 
though here, lacking any exact results for comparison, we 
cannot quantify its errors. 

Ground state energies. — Figure |l](a) shows the 
ground state condensation energies Eb = £b — (Fb\H\Fb) 
for even and odd grains (6 = and 1, respectively), 
which is measured relative to the energy of the respec- 
tive uncorrelated Fermi sea (\F ) = IIj<n c ]+ c ]- 1 Vac) 
or \Fi) = c} la+ \F }), calculated for N < 600 using both 
the canonical (E^) and g.c. {E^ 30 ) HQ approaches. The 
g.c. curves suggest a "breakdown of superconductivity" 
HQ for large d, in that E^ c = above some critical 
6-dependent level spacing df c . In contrast, the E^'s are 
(i) significantly lower than the E^'s, thus the projec- 
tion much improves the variational ansatz; and (ii) neg- 
ative for all d, which shows that the system can always 
gain energy by allowing pairing correlations, even for ar- 
bitrarily large d. As anticipated in the "breakdown 
of superconductivity" is evidently not as complete in the 
canonical as in the g.c. case. Nevertheless, some rem- 
nant of it does survive in E^ , since its behaviour, too, 
changes markedly at a 6 (and A) dependent characteristic 
level spacing 

d b (< d b C )- it; marks the end of bulk BCS " 
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like behavior for d < df , where E^ is extensive 
and the start of a f.d. plateau for d > d£ , where E^ 
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practically intensive (almost d independent) jlj]. The 
standard heuristic interpretation |L5| of the bulk BCS 
limit — A 2 /(2g?) (which is indeed reached by E^ for 
d — > 0) hinges on the scale A: the number of levels 
strongly affected by pairing is roughly A jd (those within 
A of £p), with an average energy gain per level of — A/2. 
To analogously interpret the d independence of E^ in 
the f.d. regime, we argue that the scale A loses its sig- 
nificance - fluctuations affect all uq = LUo/d unblocked 
levels within ior> of ef (this is made more precise be- 
low), and the energy gain per level is proportional to a 
renormalized coupling —Xd (corresponding to the l/N 
correction of to the g.c. BCS result). The inset of 
Fig. 0(a) shows the crossover to be quite non-trivial, be- 
ing surprisingly abrupt for E^ . 

Parity Effect. — Whereas the ground state energies 
are not observable by themselves, the parity-dependent 
spectral gaps Q.q — £2 — £q and fii = £3 — £\ are measur- 
able in RBT's experiments by applying a magnetic field 
||. Figure |l|(b) shows the canonical (rif ) and g.c. (£l£ c ) 
results for the spectral gaps. The main features of the 
g.c. predictions are |§j: (i) the spectral gaps have a mini- 
mum, which (ii) is at a smaller d in the odd than the even 
case, and (iii) Qi < Qq for small d, which was argued 
to constitute an observable parity effect. Remarkably, 
the canonical calculation reproduces all of these qualita- 
tive features, including the parity effect, differing from 
the g.c. case only in quantitative details: the minima are 
found at smaller d, and fl^ < Slf for large d. The latter 
is due to fluctuations, neglected in £ff , which are less 
effective in lowering £^ the more levels are blocked, so 
that \£fr — £^1 decreases with b. 

Wave functions. — Next we analyze the variationally- 
determined wave-functions. Each |6) can be character- 
ized by a set of correlators 

Cf (d) = {c] + c J+ c]_c 3 _) - {c) + c J+ ){c)_c^), (4) 

which measure the amplitude enhancement for finding 
a pair instead of two uncorrelated electrons in \j, ±). 
For any blocked single-particle level and for all j of an 
uncorrelated state one has Cj — 0. For the g.c. BCS 
case Cj = UjVj and the Cj's have a characteristic peak 
of width ~ A around Ef, see Fig. ||(a), implying that 
pairing correlations are "localized in energy space" . For 
the BCS regime d < A, the canonical method produces 
Cj's virtually identical to the g.c. case, vividly illustrat- 
ing why the g.c. BCS approximation is so successful: not 
performing the canonical projection hardly affects the pa- 
rameters Vj if d <C A, but tremendously simplifies their 
calculation (since the 2hq equations in (0) then decou- 
ple). However, in the f.d. regime d > d^ the character 
of the wave-function changes: weight is shifted into the 
tails far from ef at the expense of the vicinity of the 
Fermi energy. Thus pairing correlations become delocal- 
ized in energy space (as also found in pi), so that refer- 



ring to them as mere "fluctuations" is quite appropriate. 
Fig. ||(b) quantifies this derealization: Cj decreases as 
(A\£j — e f\ + B)^ 1 far from the Fermi surface, with d- 
dependent coefficients A and B; for the g.c. d = case, 
A = 2 and B = 0; with increasing d, A decreases and B 
increases, implying smaller C/s close to Ef but a slower 
fall- off far from ef- In the extreme case d 3> df pair- 
mixing is roughly equal for all interacting levels. 

To quantify how the total amount of pairing correla- 
tions, summed over all states j, depends on d, Fig. |l|(c) 
shows the pairing parameter A{,(d) = XdY) j Cj proposed 
by Ralph |||l6|, calculated with the canonical (A^) and 
grand-canonical (A^ ) approaches. By construction, 
both Af c and Af reduce to the bulk BCS order pa- 
rameter A as d — > 0, when Cj — > UjVj. Af c decreases 
with increasing d and drops to zero at the same critical 
value df c at which the energy E^ c vanishes ||, reflect- 
ing again the g.c. "breakdown of superconductivity" . In 
contrast, Af is non-monotonic and never reaches zero; 
even the slopes of AJ? and A^ c differ as d -> 0, §,§, 
illustrating that the l/N corrections neglected in the g.c. 
approach can significantly change the asymptotic d — * 
behavior (this evidently also occurs in Fig. f.b). Nev- 
ertheless, Ap c does show a clear remnant of the g.c. 
breakdown, by decreasing quite abruptly at the same d^ 
at which the plateau in E^ sets in. For the odd case 
this decrease is surprisingly abrupt, but is found to be 
smeared out for larger A. We speculate that the abrupt- 
ness is inversely related to the amount of fluctuations, 
which are reduced in the odd case by the blocking of the 
level at ef, but increased by larger A. Af increases for 
large d, because of the factor Ad in its definition, com- 
bined with the fact that (unlike in the g.c. case) the Cj 
remain non-zero due to fluctuations. 

Our quantitative analysis of the derealization of pair- 
ing correlations is complimentary to but consistent with 
that of MFF @. Despite being limited to n < 12, MFF 
also managed to partially probe the crossover regime 
from the f.d. side via an ingenious rescaling of parame- 
ters, increasing A at fixed u>d and d, thus decreasing d/ A; 
however, the total number of levels 2u>d/d stays fixed in 
the process, thus this way of reducing the effective level 
spacing, apart from being (purposefully) unphysical, can 
only yield indirect and incomplete information about the 
crossover, since it captures only the influence of the levels 
closest to ep- Our method captures the crossover fully 
without any such rescalings. 

Matveev-Larkin's parity parameter. — ML || have in- 
troduced a parity parameter, defined to be the difference 
between the ground state energy of an odd state and the 
mean energy of the neighboring even states with one elec- 
tron added and one removed: A M l = £\ -\{£q &a J r£^ m )- 
Figure |^ shows the canonical and g.c. results for Aml, 
and also the large-c? approximation given by ML, Aml = 
d/(21og(o;d/A)), where the constant a (needed because 
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ML's analysis holds only with logarithmic accuracy) was 
used as fitting parameter (with a = 1.35). As for the 
spectral gaps, the canonical and g.c. results are quali- 
tatively similar, though the latter of course misses the 
fluctuation-induced logarithmic corrections for d > d c '. 

In summary, the crossover from the bulk to the f.d. 
regime can be captured in full using a fixed- N projected 
BCS ansatz. With increasing d, the pairing correlations 
change from being strong and localized within A of Sf, 
to being mere weak, energetically delocalized "fluctua- 
tions"; this causes the condensation energy to change 
quite abruptly, at a characteristic spacing d c <x A, from 
being extensive to intensive (modulo small corrections). 
Thus, the qualitative difference between "superconduc- 
tivity" for d < d c , and "fluctuations" for d > d c , is that 
for the former but not the latter, adding more particles 
gives a different condensation energy; for superconduc- 
tivity, as Anderson put it, "more is different". 
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FIG. 1. (a) The ground state correlation energies Eb, (b) 
the spectral gaps Qb ~ £b+2 — £b and (c) the pairing pa- 
rameters Ab, for even and odd systems (b = 0, 1), calculated 
canonically (C) and grand-canonically (GC) as functions of 
d/A — 2sinh(l/A)/A r . The inset shows a blow-up of the re- 
gion around the characteristic level spacings = 0.5A and 
di — 0.25A (indicated by vertical lines in all subfigures). The 
db (a) mark a change in behavior of from ~ l/d to be- 
ing almost d independent, and roughly coincide with (b) the 
minima in , and (c) the position of the abrupt drops in Aj,. 
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tive pairing parameter proposed in Bj, turns out to be 
identically equal to A& for the ansatz \b). 




FIG. 2. The pairing amplitudes Cj of Eq. QJ). (a) The 
dashed line shows the g.c. BCS result; pair correlations are 
localized within A of ef- Lines with symbols show the canoni- 
cal results for several d; for d < d$ « 0.5 A, the wave functions 
are similar to the BCS ground state, while for d < d$ weight 
is shifted away from ef into the tails, (b) For all d, GJ 
shows linear behavior far from ef- For larger d the influence 
of levels far from ef increases. 




FIG. 3. The canonical (solid), g.c. (dashed) and analytical 
(dotted line) results for the parity parameter A M l |j ■ 
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